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We calculated the superconductivity properties of alkali metals under high pressure using the 
results of band theory and the rigid-mumn-tin theory of Gaspari and Gyorffy. Our results suggest 
that at high pressures Lithium, Potassium, Rubidium and Cesium would be superconductors with 
transition temperatures approaching 5 — 20K. Our calculations also suggest that Sodium would 
not be a superconductor under high pressure even if compressed to less than half of its equilibrium 
volume. We found that the compression of the lattice strengthens the electron-phonon coupling 
through a delicately balanced increase of both the electronic and phononic components of this 
coupling. This increase of the electron-phonon coupling in Li is due to an enhancement of the s-p 
channel of the interaction, while in the heavier elements the p-d channel is the dominant component. 



I. INTRODUCTION 

Neaton and Ashcrofti predicted that at high pres- 
sures Li forms a paired ground state. Subsequently, 
Christensen and Novikov^ showed that fee Li un- 
der increased pressure may reach a superconducting 
transition temperature T c =50— 70K. This prediction 
has been supported by Shimizu et a& and Struzhkin 
et al* who reported measurements of superconduc- 
tivity in compressed Li with a T c ranging from 9K 
to 20K. The overestimate of T c reported in Ref[2] 
is probably due to an error by these authors as we 
discuss in Sec. III. Using the rigid-muffin-tin approxi- 
mation(RMTA) formulated by Gaspari and Gyorffy D , 
we performed calculations for two other alkali met- 
als, K and Rb, and predicted that at high pressures 
they both would be superconductors with transition 
temperatures approaching 1QK 6 . Before these works, 
the only known superconducting alkali metal was Cs, 
which becomes superconducting above 7.0 GPa with 
a transition temperature of about 1.5K 7 Similar the- 
oretical results for Li, K and Cs were also reported by 
Ashcroft 8 , Profeta et ali£, Kasinathan et aliiS, Tse et 
alii and Stocks et a\J£ 

In this paper, we have extended our study to Li, 
Na and Cs using again the RMTA. Our calculations 
demonstrate that Li and Cs display superconductiv- 
ity at high pressure above 15 GPa and 3.5 GPa re- 
spectively. Our results also showed the lack of su- 
perconductivity for Na up to 90 GPa. We compare 
these new calculations with our previous results of K 
and Rb to show a complete picture of superconduc- 
tivity properties of alkali metals under high pressure. 
Our calculations indicate that the s-p channel of con- 
tribution to the Hopheld parameter rj dominates Li 
under high pressure, while the p-d channel contribu- 
tion is the major reason that K, Rb and Cs become 
superconductors under high pressure. 

The massive structural phase transitions of all 
alkali metals have be extensively investigated dur- 
ing the past few decades, n i 13 i 14 i 15 i 16 i 17 i 18 i 19 All al- 
kali metals were found to be stable with bee struc- 



ture under ambient pressure and will transform to the 
fee structure at about 20, 65, 11.5, 7, and 2.3 GPa 
pressure for Li, Na, K, Rb and Cs respectively. In 
this paper, our investigations are focused on the bee 
and fee structures within the above pressure ranges 
which correspond to volume changes as large as 40% 
from equilibrium. The conclusions we draw in this 
paper are valid for the structure which is stable in 
the experimental range of pressure. 

Tomita et alSI recently reported experimental re- 
sults for Li, Na and K. They confirmed the super- 
conductivity of Li above 20 GPa at temperatures 
reaching 15 K. They also pointed out the absence 
of superconductivity in Na and K for pressure up 
to 65 and 43.5 GPa respectively without specifying 
whether the measurements were extended to the fee 
phase. However, for K, our calculations and those 
of Profeta et ai£ show superconductivity in the fee 
phase. Further experimental work may be required. 



II. THEORY AND COMPUTATIONAL 
DETAILS 

McMillan's strong coupling theorjiSi defines an 
electron-phonon coupling constant by: 



A = 



M < uj 2 > 



N(e F ) <I 2 > 
M < uj 2 > 
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where M is the atomic mass; r\ is the Hopfield param- 
ete»2£; N(ep) is the total density of states(DOS) per 
spin at the Fermi level, ep\ <I 2 > is the square of the 
electron- ion matrix element at ep; and <u>> is the 
average phonon frequency. The <I 2 > is determined 
by the RMTA approximation given by Gaspari and 
Gyorffj£(GG) formula: 
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where Si are scattering phase shifts, Ni(ep) is the Ith 

component of the DOS per spin, and is the free 
scatterer DOS. 

There are two important papers that address the 
question of the spherical approximation in the GG 
theory. The first one by John2& shows that the GG 
formula is exact for cubic systems with one atom in 
the unit cell and with l+l up to 2. This validates our 
results for Li and Na where the p-d and d-f chan- 
nels have negligible contributions. In addition, But- 
ler et ah 24 have obtained a generalization of the GG 
formula for l>2 that is exact for systems in which 
all atoms sit at sites having cubic symmetry. Their 
expression contains certain cross terms which they 
showed have small or canceling contributions to the 
value of 77. Therefore, our calculations for K, Rb, and 
Cs which have strong diagonal contributions from p-d 
but not d-f scattering are also reliable. 

In order to determine the quantities entering 
Eq. [21 we performed Augmented Plane Wave ( AP W) 
calculations of the band structures and total ener- 
gies of the targeted alkali metals in the local density 
approximation(LDA) following the Hedin-Lundqvist 
prescription 2 ^ in both bcc and fee structures, for a 
wide range of pressures. Since the LDA underesti- 
mates the equilibrium lattice constant of alkali met- 
als by approximately 3% we applied a uniform shift 
to our results for E(V) so that at the experimental 
equilibrium volume we have the minimum energy. To 
accomplish this the pressure shifts are 1.5, 0.4, 0.4, 
0.4 and 0.25 GPa for Li, Na, K, Rb and Cs respec- 
tively. We also used the APW results to determine 
the pressure variation of the bulk moduli B. The B 
values were used to compute the average phonon fre- 
quency which has been taken to be proportional to 
the product of bulk modulus and volume 26 : 



< lo 2 >= CB{V)V 



1/3 



(3) 



where the constant C is determined by the <uj 2 > and 
the Debye temperature Ou under normal pressure 
using the formula: 



= lei 



(4) 



The approximation made in using Eq.[3]and Eq.0] 

gives reasonable values of A and leads us to believe 
that our results should differ only quantitative from 
direct phonon evaluations via linear response theo- 
ries. The transition temperature for superconductiv- 
ity is given by the McMillan equation^!: 



t <UJ> r 
T c = 1 2 exp[ 



-1.04(1 + A) 
A-/i*(l + 0.62A) J 



(5) 



III. RESULTS 

Before we present our detailed results, we wish 
to bring to the attention of the reader that in our 
opinion there must be an error in the evaluation of 
77 for Li by Christensen and NovikoviS The error is 
that in Eq.^they must have used N(cf) for two spins 
instead of one as is the correct implementation of the 
Mcmillan theory. Fig. ^left panel) demonstrates a 
clear discrepancy between our results and those of 
Ref 2. Dividing the 77 of Ref 2 by two gives a perfect 
agreement with our results as shown in Fig. fright 
panel) and as we discuss below lowers the value of A 
and brings T c close to the measured value. 



Pressure(GPa) 

6.5 2.9 1.3 



Pressure(GPa) 

6.5 2.9 1.3 




where /j* is the Coulomb pseudopotential, that we 
have set equal to 0.13. 



FIG. 1: Comparison of 77 between our results and Chris- 
tensen and Novikov for fee Li in the left panel, and the 
corrected result in the right panel. The vertical dashed 
line indicates the transition from bcc to fee. 

To demonstrate the effect of high pressures on the 
density of states at the Fermi level for alkali metals 
we plotted the total density of states and its angu- 
lar momentum decompositions for fee Li and Cs un- 
der ambient and high pressures in Fig. [2] For Li, we 
note that both N(£f) and its angular momentum de- 
compositions decrease under pressure. On the other 
hand, for Cs, we observe a decrease of the s and p- 
like density of states with a remarkable increase of 
the d-like density of states N^ep) under high pres- 
sure. This increase of the d density of states at ep 
makes the largest contribution to the large value of 
77 at large pressure as discussed below. We also ob- 
served the same trend as in Cs for the alkali metals K 
and Rb. In contrast, for both the bcc and fee phases 
of Li and Na, the density of states decreases as the 
lattice is compressed. This fact gives a partial ex- 
planation of the absence of superconductivity under 
high pressure for Na. For Li, however, the increase in 
the matrix elements <I 2 > is strong enough to ensure 
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the overall large increase in 77, which we will discuss 
later. 
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FIG. 2: Total density of states and angular momentum 
decompositions of fee Li and Cs under normal and high 
pressure. The solid and dashed lines denote the high and 
ambient pressure respectively. 

We show the ratios Ni(eF)/N(ep) as a function 
of pressure in Fig. [3] where the vertical line indicates 
the pressure where the transition from bcc to fee oc- 
curs in the experiments. These ratios are crucial in 
the determination of 77. It is important to note that 
the ratio Nd(eF)/N(eF) of the heavier alkali metals 
K, Rb and Cs increases rapidly as a function of in- 
creasing pressure. The buildup of the d-like DOS 
under high pressure causes the large values of 77 for 
the heavier elements under high pressure shown in 
Fig. 0] For Li and Na, the d ratio is very small and 
the s ratio dominates. We note that for Na, the p 
ratio decrease with increasing pressure. 

In Fig. we show the total 77 and its contribu- 
tions from s-p, p-d and d-f scattering channels as 
a function of pressure in both the bcc and fee struc- 
tures. 77 in all five alkali metals in the bcc and fee lat- 
tices increases significantly with increasing pressure. 
Among all alkali metals, Li has the second largest 
increase of 77 with increasing pressure, which is one 
of the reasons why Li is a superconductor under high 
pressure. Fig. 0] shows that different elements have 
different scattering channels as the dominant contrib- 
utors to the total 77. More specifically for Li and 
Na, the 77 contribution of the s-p channel increases 
rapidly with increasing pressure while the other two 
contributions stay around zero. For K and Rb, when 
pressure is increased the 77 contributions from both 
the s-p and p-d channels increase quickly. For Cs, the 
p-d channel contributes the major increase of 77 with 
pressure increasing. For all elements, the 77 contribu- 
tion of d-f channel is always very close to zero. The 
graph also demonstrates that the largest portion of 



005 FIG. 3: Angular momentum decomposed DOS divided 
by the total DOS at ef, the filled square, circle and dia- 
mond symbols denote s, p and d states respectively. 



the 77 increase in Li and Na under high pressure is 
contributed by the s-p channel. Na shows the largest 
value of 77 among all the alkali metals when pressure 
reaches 90 GPa. However, as we discuss below Na has 
a very large value of <u>> which lowers significantly 
the value of A. 
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FIG. 4: Total r\ and its contributions from the s-p, p-d 
and d-f channels, the filled circle, square, diamond and 
triangle symbols denote total, s-p, p-d and d-f chanels 
respectively. 

Fig. [5] shows the increase of the bulk modulus B 
as a function of pressure. We note that the bulk 
moduli of Li and Na increase dramatically with in- 
creasing pressure. The bulk modulus increase of K, 
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FIG. 5: The bulk moduli B as a function of pressure, 
the filled circle and square symbols denote bcc and fee 
respectively. 



FIG. 6: Average phonon frequency < uj > as a function 
of pressure. The filled circle and square symbols denote 
bcc and fee respectively. 



Rb and Cs is around 1/3 to 1/10 of that for Li and 
Na. Fig. |B| shows the prefactor, <ui>, in the McMil- 
lan equation. We note that the <uj> for Li and Na, 
in both bcc and fee structures, increases much faster 
than in the other alkali metals. This is the other 
reason that Li has larger T c than K, Rb, and Cs de- 
spite the fact that A of K, Rb, and Cs can be larger 
than that of Li under pressure. Combining the val- 
ues of 77 and M<uj 2 >, we determined the electron- 
phonon coupling constant A using Eq. ^ In Fig. [7J 
we show A as a function of pressure. It is evident 
that for the alkali metals Li, K, Rb, and Cs, under 
pressure, A reaches large values suggesting that these 
metals could display superconductivity. To quantify 
our predictions for superconductivity, we have calcu- 
lated T c for all five elements in both the bcc and fee 
structures at high pressure using the McMillan equa- 
tion with Coulomb pseudopotcntial values jj,* = 0.13 
and 0.20. Fig.[H]shows the T c of the five alkali metals 
in the fee structure as a function of pressure. As one 
might expect from the large values of the electron- 
phonon interaction A, the elements Li, K, Rb, and 
Cs are predicted to be superconductors with T c larger 
than 5 K for pressures higher than 20, 11, 7, and 3.5 
GPa respectively. Clearly our results for the value 
of T c are sensitive to the value of /z* and increase of 
fi* to 0.2 suppresses T c by a factor of four for Li and 
about 20% to 40% for K, Rb and Cs, as shown in 
Fig.0 

Our calculations suggest that Na does not display 
superconductivity because the electron-phonon cou- 
pling constant A remains small (A « 0.25) for com- 
pressions up to 90 GPa as shown in Fig. [7\ The 
absence of superconductivity for Na is determined 
by delicate balance of the increasing electronic and 
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FIG. 7: Electron-phonon coupling constant A as a func- 
tion of pressure. The filled circle and square symbols 
denote bcc and fee respectively. 



phononic components of the electron-phonon cou- 
pling constant. Na has a r\ which increases faster 
than in K, Rb and Cs, but slower than in Li. How- 
ever, the increase of r\ in Na is mostly canceled out by 
its M<u) 2 > which increases much more rapidly than 
in K, Rb and Cs. Therefore, Na does not display 
strong enough electron-phonon coupling to become a 
superconductor under high pressure up to 90 GPa in 
both the bcc and fee structures. 
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FIG. 8: Transition temperature T c as a function of pres- 
sure, the filled circle and square symbols denote bee and 
fee respectively. 



IV. CONCLUSIONS 

As it can be seen from Eq.^ the increase of A may 
be caused by either an increase of r) or a decrease of 
M<oj 2 >, or a combination of these two factors. For 
K, Rb, and Cs, because the increase of the bulk mod- 
ulus under pressure is slow, the increase of 77, which 
is mainly caused by an increase of N(ep), dominates 
the determination of A. Therefore, K, Rb, and Cs 
have a large A under high pressure, which results in 



the prediction of superconductivity. The increase of 
A in Li under pressure is due to the increase of <I 2 > 
rather than N(ep), which is actually decreasing un- 
der pressure, and overcomes the increase of M<u> 2 > 
which in fact helps by increasing the prefactor <c<j> 
in the T c equation. Therefore, Li would be a super- 
conductor under high pressure. Na is different from 
Li and from the other alkali metals because its large 
r\ is canceled by the increasing bulk modulus under 
pressure which reduces the value of A. Therefore, a 
delicate balance of rj and M<uJ 2 > results in Li being 
a superconductor but not Na. Our calculated T c for 
Li showing in Fig [5] is in very good agreement with 
experiment £ 

It should be mentioned here that Tse et alii ar- 
gued that the softening of the TA phonon branch 
near the X point causes the occurrence of supercon- 
ductivity in Li. Our view is that this is an over- 
simplification and is not a sufficiently quantitative 
explanation. On the other hand the RMTA takes 
an average in all directions on the Fermi surface and 
not just in one high symmetry point. Regarding the 
mechanism of superconductivity in K, Rb, and Cs 
we suggest that this is due to the increased d-like 
character of the wave functions at ep at high pres- 
sures. This enhanced d character near has also 
been documented by other authors^ and by experi- 
ment. We used it here to justify the RMTA which 
is successful in transition metals as shown by Papa- 
constantopoulos et aliSS 
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